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. Abstract. BRAHMS has the ability to study relativistic heavy ion coUisions from 

I the final freeze-out of hadrons all the way back to the initial wave-function of the gold 

' nuclei. This is accomplished by studying hadrons with a very wide range of momenta 

Oh, and angles. In doing so we can scan various phases of QCD, from a hadron gas, to a 

quark gluon plasma and perhaps to a color glass condensate. 
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1. Introduction 
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^ ■ The purpose of RHIC is to map the phase structure of QCD. So far the community has 
O ' concentrated on AuAu, dAu and pp collisions at ^JsNN = 200 GeV in the hope of finding 

O 



the quark gluon plasma. BRAHMS' special contribution has been to study the hadrons 
produced in these collisions over a broad range of pr and rapidity p. A great deal of 
^ ■ evidence supporting the creation of partonic matter in AuAu collisions was presented 
Q ' at this conference. However QCD is a rich theory that probably has many phases. It 
has been suggested that when viewed by a fast probe a heavy nucleus may resemble a 
sheet of highly correlated gluons called the Color Glass Condensate |2]. After a AuAu 
collision the fields generated by the color charges on the two sheets of gluons would break 
up into partons which one would expect to approach chemical and kinetic equilibrium 
while rapidly expanding in both the longitudinal and transverse directions. Eventually 
the partons must hadronize and after further rescattering the hadrons freeze-out. 

This scenario is speculative and our evidence is both incomplete and somewhat 
indirect. Nevertheless in this paper we will attempt to map out this evolution by 
starting from the final state and working our way backwards. The RHIC experiments 
have a beautiful complimentarity but we will report only on BRAHMS' data, with an 
emphasis on recent results. In particular we will discuss: 

• Multiplicity distributions from dAu collisions jS] which reflect entropy production. 

• Particle spectra from AuAu which give information regarding; 

— kinetic freeze-out via blast wave fits to px spectra, 

— chemical freeze-out via fits to particle ratios, 

— initial pressure and longitudinal flow from pion dN/dy distributions. 

• High pt suppression j5] which is sensitive to the early density of color charges. 

• The ratio of dAu and pp spectra at high rapidities jB] which give information on 
the Au wavefunction. 
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2. Global Observables 

Multiplicity distributions are sensitive to all stages of the collision and can be used 
to measure the total production of entropy. Figure [T] shows our dN/drj results for 
minimum-bias and central dAu collisions 3J (for AuAu see [3 Ej)- Panel (c) shows 
the ratio of the 0-30% and 30-60% samples normalized by the number of participants. 
The ratios appropriate for Au- and d-participant only scaling are indicated by the left 
and right arrows. Particle production away from mid-rapidity appears to follow the 
participant scaling of the respective fragment. In the deuteron fragmentation region we 
see very similar yields to lower energy data. This phenomenon is known as "limiting 
fragmentation" P El CD]- The HIJING and AMPT models are close to the data 
[TT| 1121 ESI- Note that the saturation model results have been updated since the 
conference with a better centrality determination and an increase of the saturation 
scale from Ql = 0.25 to O.^AGeV"^ These new calculations are close to the data. 




^ n n 

Figure 1. dN/drj distributions from dAu collisions 3 for a) minimum- bias and 0-30% 
central events and b) 30-60% central events, c) Scaled multiplicity/participant ratio 
R. The left (right) arrows show corresponding values for Au- (d-) participant scaling. 



3. Particle Spectra 

The distribution of particles in rapidity and px may give information on the transverse 
and longitudinal flow while the mix of different kinds of particles may tell us about 
the "quark chemistry" of the system. Our AuAu spectra are summarized in Fig. 
which shows the rapidity densities, dN/dy, and the mean transverse momenta, (pr) , 
for TT, p and p as a function of rapidity. Both quantities are estimated using fits to 
the spectra in narrow regions of rapidity |4j. For vr. A; and p the yields peak at y=0 
and drop significantly at higher rapidities. The 7r+ and tx" yields are nearly equal while 
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Figure 2. Rapidity densities (top) and mean transverse momentum (bottom) as a 
function of rapidity. Errors are statistical. The yields of all the produced particles are 
well described by Gaussian fits, shown by the dotted lines. The beam rapidiity is 5.4 



an excess of over K~ is observed that increases with rapidity. The lower panel of 
Fig. 121 shows the rapidity dependence of (pt) ■ There is no significant difference between 
particles and their antiparticles. In general, the rapidity dependence of {px) increases 
with mass suggesting that transverse flow drops with increasing rapidity. Using the net 
proton yield combined with baryon conservation and some assumptions on the neutron 
and hyperon yields allows us to estimate the total energy liberated by the stopping of 
the baryons. We find AE = J EdN = 25 ± 1 TeV, or 75 GeV per participant pilT^. 

3.1. Rapidity Dependence of Kinetic and Chemical Freeze-out 

Is there one source or many in high energy heavy ion collisions? We have investigated 
this question by fitting our spectra and particle yields at several different rapidities to 
blast wave and chemical models ^Hl Ej- At y=0 we see a very slow change of the 
freeze-out parameters with centrality so we shall consider only central data here. The 
left panel of Fig. El shows the regions of temperature T and transverse velocity of the 
surface [3s that are consistent with our data sets at y=0,l,2 and 3. As the rapidity 
increases j3s decreases while T increases. This may be because the equation of state of 
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the matter is changing with rapidity. If the number of degrees of freedom decreases one 
would expect the temperature to increase. The right panel of Fig. El shows the results 
of a chemical analysis versus rapidity. As y increases both the baryo-chemical potential 
and (to a lesser extent) the chemical freeze-out temperature increase. Again this may 
suggest that the system has fewer degrees of freedom at higher rapidities. 




Figure 3. Preliminary blast wave (left) and thermal (right) fits to our data at various 
rapidities. The curves are the one sigma contours. Note that as expected the thermal 
temperatures are always higher than the kinetic freeze-out temperatures. 



3.2. Bjorken and/ or Landau Hydrodynamics 

We now turn to the longitudinal flow which may be sensitive to the initial pressure in 
the system and possibly the equation of state. Bjorken proposed ^Hj that away from 
the fragmentation regions, y ~ ±4 at RHIC ^3], the system produced by heavy ion 
collisions should be boost invariant, i.e. independent of rapidity. This assumption is 
pervasive in the theoretical literature. Such an expansion is the fastest possible one 
in the longitudinal direction. If the expansion is slower than the Bjorken limit then 
freeze-out occurs later and it may be easier to explain the fact that the pion HBT radii 
are the same in the "sidewards" and "outwards" directions. 

For \y\ < 1 all of our data are consistent with Bjorken's proposal. However 
looking globally at Fig. |21 there is a clear breakdown of boost invariance. This is 
most noticeable in the particle yields but it is also true that the (pr) of the kaons 
and antiprotons falls significantly with rapidity. Clearly a full understanding of the 
longitudinal dynamics would explain the vr, k and p data. However because the pions 
dominate both the multiplicity and transverse energy, Et = \I'Pt + distributions 
focussing on the pions is a good start to a complete description of the longitudinal 
flow. Landau developed an analytic model of relativistic hydrodynamics undergoing an 



Scanning the Phases of QCD 



5 



isentropic (constant entropy) expansion governed by an equation of state . This idea 
was extended by Carruthers et al pi? to pion rapidity distibutions by assuming that 
pions are massless and that their pt and rapidity distributions approximately factorize. 
Under these conditions dN/dy is a gaussian with a width given by 

where is the nucleon mass. 

This model was able to give a reasonable description of the pion distributions from 
pp collisions at various energies. The assumptions of the model are not entirely met for 
our data since = 0.3- (pt) at y=0 and {px) drops by 10% from y=0 to y=3. Another 
difference between our data and the Landau assumptions is that we do not observe full 
stopping. Nevertheless the agreement between this very simple model and our data is 
rather good. FigurelUJa) shows dN/dyiji) and Landau's prediction for ^^snn = 200 GeV 
using Eq. [T]with the condition that the integrals of these Gaussians must be equal to 
the full-space yields estimated from the data. A discrepancy of ~ 5% is observed 
{dLandau = 2.16). Figure El^b) shows a compilation on pion widths from AGS to RHIC, 
The difference between theory and data is at most 10%. The logarithmic growth of 
the rapidity width with ^/sNN is in contrast to the linear increase of the multiplicity 
with y/SNN [21]. It is all the more striking considering that the degree of transparency 
drastically changes from AGS to RHIC energies |15j . 



4. High Pt suppression and energy loss 

The most exciting heavy ion news of 2003 was the discovery that high px suppression 
in AuAu collisions is a result of the hot and dense medium produced in these collisions 
rather than a depletion of hard partons in the Au nucleus itself [221 ■ We quantify this 
effect by normalizing our spectra to pp data using the nuclear modification factor; 

Here (Ncoii) is the average number of nucleon- nucleon collisions in each event. BRAHMS 
has the unique ability to study this effect over a wide rapidity range 0. Figure El shows 
Raa ioT TT^ at y=2.2 for AuAu and dAu collisions. For central dAu collisions there 
is already some suppression at y=2.2 although it is not as strong as in central AuAu 
coUisons [221 IS] • Note however that this measurement relies on an extrapolation for the 
pp reference based on PHENIX measurements, see 



5. The Initial Gold Wavefunction 



Finally it has been suggested that when viewed by a fast probe a heavy nucleus may 
form a new phase of QCD, the Color Glass Condensate Pj. In such a state, two soft 
gluons will often fuse together to produce one harder one. This effect results in a low 
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Figure 4. (a) Pion dN/dy and Landau's Gaussian prediction at ^Jsmn — 200 GeV 
from Eq. ^ The errors are statistical, (b) <yN(-n)/^ Landau as a function of y^ijvw 
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Figure 5. Preliminary tt^ nuclear modification factors at y—2.2 for AuAu and dAu 
collisions. The dashed and dotted lines represent the value expected if the yield scales 
by (Ncoii) or (Nparticpant) respectively. The bands show our systematic errors on these 
quantities 



Pt suppression for dAu collisions compared to pp ones. BRAHMS can tune the speed 
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of the probe by looking at fragments that emerge over a wide range of rapidities. Since 
we observe fragments close to the axis of the deuteron beam we are observing the initial 
parton density in the gold nucleus. The partons momentum fraction, probed by a given 
fragment is given by x = e'^pxl ^/snn- As y increases, x decreases and if the Color 
Glass scenario is correct we should see a shift of low px particles to high pt with an 
overall reduction in yield. 

Figure El shows RdAu for minimum-bias events as a function of pt and rj. We 
use 7] instead of y because we wish to combine all particles together and so improve 
our statistical accuracy. RdAu rises with px and falls with rj. The systematic errors 
in RdAu range from < 10 — 15%. At midrapidity, RdAu goes above 1. This so-called 
Cronin enhancement [23 has been attributed to multiple scattering of the incoming 
partons during the collision. At 77 = 1 the Cronin peak is not present and at more 
forward rapidities (77 = 3.2) the data show a suppression of the hadron yields. At low 
VT) RdAu is close to the ratio of charged-particle pseudorapidity densities in dAu and 
pp collisions [SlEll- Saturation effects should increase with the thickness of nuclear 
material traversed by the incoming probe and indeed we see a greater suppression for 
more central collisions, see Fig. El and refs. [23 IE] • 
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Figure 6. Nuclear modification ratio for charged hadrons versus pr and rj. Systematic 
errrors are shown with shaded boxes. The band around unity shows the error on (Ncou)- 
Dashed Hnes from px — — 1.5 GeV/c show the ratio (jv^ '^dli'jdi^ipp) ■ 



6. Summary and Conclusions 

For dAu collisions we see a significant asymmetry in dN/drj with a peak at rj = —2 
(i.e. on the Au side of the collision) and a slight shoulder at 77 = +2. This suggests 
significant rescattering within the dAu system since these peaks are far away from the 
Au and d beam rapidities. In the fragmentation regions the multiplicity scales with the 
number of (the Au or deuteron) particpants. Using our spectra of identified particles we 
have found that the rapidity distributions of all the produced charged particles in AuAu 
collisions are Gaussian. There is no large rapidity plateau but the data are consistent 
with boost invariance for \y\ < 1. The width of our pion distribution (and a large range 
of lower energy data) is consistent with Landau's picture of isentropic fluid dynamics. 
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Blast wave analyses of our data show a decrease in the surface velocity [3s and 
an increase in the kinetic freeeze-out temperature with increasing rapidity. Similarly 
chemical analysis of our particle yields hint that both the baryo-chemical potential and 
the chemical freeze-out temperature increase with rapidity. One could interpret this in 
terms of the system becoming less partonic (with consequently fewer degrees of freedom) 
at higher rapidities. The thermal analysis allows us to make a rough estimate of the 
total energy in the produced particles. This yields 25 ± 5 TeV compared to 25 ± 1 TeV 
computed from integrating the energy in the net protons [TH] . 

We see evidence for jet quenching in AuAu collisions at both y=0 and y=2.2. 
The simplest explanation for these data is that fast partons lose energy traversing a 
deconfined system of quarks and gluons In the forward region we see a low 

suppression of the yield from dAu collisions which increases with r] and centrality. This 
may be a result of a saturation in the yield of low momentum gluons in the gold nucleus 
and hints that the Color Glass Condensate may represent the high energy limit of QCD. 
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